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Abstract 

A new theoretical model is introduced to describe heterodyne mixer conversion 
efficiency and noise (from thermal fluctuation effects) in diffusion-cooled superconducting 
hot-electron bolometers. The model takes into account the non-uniform internal electron 
temperature distribution generated by Wiedemann-Franz heat conduction, and accepts for 
input an arbitrary (analytical or experimental ) superconducting resistance-versus- 
temperature curve. A non-linear large-signal solution is solved iteratively to calculate the 
temperature distribution, and a linear frequency-domain small-signal formulation is used to 
calculate conversion efficiency and noise. In the small-signal solution the device is 
discretized into segments, and matrix algebra is used to relate the heating modulation in the 
segments to temperature and resistance modulations. Matrix expressions are derived that 
allow single-sideband mixer conversion efficiency and coupled noise power to be directly 
calculated. The model accounts for self-heating and electrothermal feedback from the 
surrounding bias circuit. 


1 . Introduction 

In the last few years, the superconducting hot-electron bolometer (HEB) has attracted 
much interest for use as the. mixing element in low-noise submillimeter wave spectrometer 
front-ends for astrophysics. Instruments are now being developed for such missions as the 
National Aeronautics and Space Administration's airborne SOFIA platform, and the 
European Space Agency's spacebome FIRST telescope. One of the key properties that 
make HEB mixers interesting is their ability to operate significantly above the 
superconducting gap frequency, which is a materials-related parameter that currently limits 
competing SIS mixers to frequencies below about 1200 GHz. Even though this is a new 
heterodyne detector that with one exception [1] has not been used for practical astrophysical 
observations (primarily because frequencies above 1 THz cannot be observed from ground- 
based observatories), rapid improvements in measured receiver noise temperatures [2-6], 
and intermediate frequency (if) bandwidths [6,7], have been reported, primarily in 
anticipation of upcoming missions such as those mentioned above. Despite these recent 
encouraging results, further improvements in the noise, and indeed the overall optimization 
of this mixer, will require an accurate and detailed theoretical model which accounts for the 
unique operation of this device. 

Superconducting HEB mixers have mostly been analyzed using lumped element models 
that were originally developed for semiconductor bolometers [8-12], In this approach the 
bolometer is represented as a single (i.e. “lumped”) heat capacitance that is connected to a 
thermal bath via a heat conductance. While this theoretical model has significant merit, 
especially for phonon-cooled devices that have a uniform electron temperature [13], it does 



not accurately take into account the non-uniform temperature distribution in the diffusion- 
cooled case [14]. The most obvious complication lies in how to determine an appropriate 
value for the resistive transition dR/dT, which enters as a square in the lumped-element 
expression for the mixer conversion efficiency. Inclusion of a non-uniform temperature 
complicates the model enough, that a closed-form analytic solution for mixer conversion 
and noise cannot be found. Several methods to avoid this problem have been previously 
investigated. In one approach, a reduced or "effective" value for dR/dT is calculated 
numerically taking a non-uniform temperature into account, and this value is then used in 
the regular lumped-element mixer model [15]. Another solution is the so-called "hot-spot" 
model, where the device resistance-versus-temperature (RT-) curve is approximated by a 
step function [16-17]. A third method (previously employed by us) is to discretize the 
bolometer into segments, and do a brute-force numerical time-domain mixing simulation, 
calculating heat dissipation and heat flow between the segments within each time step. 
Results of this method are compared to experiments in [18]. The new frequency-domain 
method described in this paper has advantages compared to all of these previous efforts: 1) 
It allows mixer conversion efficiency and output noise to be calculated without making 
simplifications in the RT-curve; 2) The mixer performance can be determined by 
straightforward matrix multiplications and inversions, without any equivalent of the time- 
step dependent instabilities that can occur in time-domain simulations; 3) The dominant 
noise mechanism is included in the model; and 4) Self-heating and electrothermal feedback 
are fully included. There are also some issues not currently addressed this new technique. 
For example, it does not include a method for calculating rf or if saturation effects, and it 
does not include several effects that are particular to superconductors, such as the 
coherence length and the differences in heat conductance between superconductors and 
normal metals. Future inclusion of these effects would affect the first, non-linear, part of 
the theory. The small-signal part of the theory uses parameters calculated from the non- 
linear part, but may otherwise remain unaffected. 


2 The Mixer Model 

The type of HEB device we are modeling consists of a submicron-length strip 
(microbridge) of a superconductor between two normal metal leads. The pads are assumed 
to be at a bath temperature T am i,. Heat dissipation in the microbridge raises the temperature 
of the center of the device relative to the ends, creating a non-uniform temperature profile. 
The bridge length L is assumed short enough that the dominant heat conduction out of the 
bridge is by diffusion of hot electrons to the cold normal metal pads, rather than by 
electron-phonon interaction. 


2. 1 Steady State / Large Signal Model 

The large-signal temperature distribution in a diffusion-cooled HEB mixer is a result of 
the rf heating caused by the local oscillator, by the heating from the mixer dc bias current, 
and by internal heat conduction. The heat conduction is assumed to be proportional to the 
local temperature according to the Wiedemann-Franz law: 

(Eq.l) 

; o 

where O is heat flow intensity, T is the local temperature, and Go is the heat conductivity at 
a reference temperature Tq. The relation leads to the following one-dimensional integral 



equation, where the heat flow towards the ends of the bolometer at a particular distance x 
from the center is equated to the amount of power dissipated inside of that distance: 
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(Eq.2) 


where w,h,L are the width, thickness and length of the bolometer, Prfis the amount of local 

oscillator power dissipated in the bridge, 1q is the dc current, and p is the local resistivity in 
the bolometer. Since the device film resistivity is usually a non-linear function of 
temperature and location, no general closed-form solution can be found to this equation. 
However, a numerical solution can quite easily be found through a corresponding finite- 
difference formulation: 
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where the integral in the previous expression is approximated by a summation Ax S k , with 

1 1 
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1 1=2 1 

The finite-difference formulation can be seen as a sub-division of the bolometer into N 
segments* each with its own temperature Tk, as is illustrated in Fig.l . Only half of the 
bolometer needs actually be considered because of symmetry. The finite-difference 
formulation is computationally fast even when used iteratively to find numerical values of 
various parameters such as the dc current 7 0 or the temperature 7J of the device center 
segment. 
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Fig.l; Discretization used in the large- and small-signal models. Because of symmetry only 
half of the bolometer is analyzed. Therefore the total length of the device is L = 2 • N • Ar . 


2.2 The Mixer Embedding Circuit 

As shown in Fig.2, the discretized bolometer is connected to a Norton equivalent circuit 
to represent the dc bias and intermediate frequency if embedding circuits of the mixer. Rq 



and / 0 are the steady-state (large signal) resistance and bias current of the device, while 

A R and A I are the small-signal modulations representing the modulation of the dissipated 
rf power caused by the combined signal and LO. 
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Fig.2; A Norton equivalent circuit is used to represent the dc bias and intermediate 
frequency circuits. 


In the small-signal limit the embedding circuit imposes the following relationship between 
resistance and current, which must be included in die subsequent mixer theory: 
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2.3 Mixer Conversion Efficiency Model 

In bolometric mixing, the device resistance is varied by the modulation in rf power 
dissipation that is caused by the beating between signal and local oscillator (LO). This 
modulation in resistance is then detected by an external circuit. In superconducting HEB 
mixers that are used at frequencies above the superconducting gap frequency, the //heating 
can be regarded as uniform throughout the device. One can also assume that the signal and 
local oscillator frequencies are higher than the highest response frequency of the device. In 
this situation it is convenient to consider the power dissipation from the local oscillator to 
be constant with time and uniformly distributed throughout the device, and to represent the 
modulation due to the signal by a uniformly distributed heating function that is modulated at 
the intermediate frequency, if. In the small-signal limit this means that the modulation of the 
device resistance will occur only at the if 

The resistance modulation A R is the sum of the modulations A R k in all the segments: 
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(Eq.7) 




where AT k is the temperature modulation in segment k. The total heat modulation in 

segment k is the sum of the contributions from rf power dissipation and from the 
embedding circuit: 


= A Prf, k + A P embtk = A P iJk + 2%/ 0 A7 + I$AR k (Eq.8) 

The frequency-domain model illustrated in Fig.C is used to calculate the effect that the 
heating modulations have on the temperatures. In the small-signal limit the thermal model is 
linear, so heat conductances, heat capacitances and heat sources can be handled by the same 

complex-notation (i.e.y'cy-notation) that is commonly used for small-signal electrical circuits 
and components. The C k ’s in Fig.C are the electronic heat capacitances of the individual 
segments, which are proprtional to the volume of the segment and to the local temperature 
T|<. The G k s are the segment-to-segment heat conductance, which according to the 
Wiedemann-Franz law is proportional to the electrical conductance and to the local 
temperature.To handle the many device segments, a matrix notation is introduced (with 
square brackets denoting matrices and vectors). The relationship between heat dissipation 
and temperature can be written as an analogue to Ohm's law: 

[AJi]*=[Au| a -[AP / ] / (Eq.9) 


Where [A] is calculated from its inverse, which can be derived according to Fig.3 as : 
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Fig.3; Thermal circuit model, showing heat dissipation sources (APk), segment heat 
capacitances (Ck) and segment-to-segment heat conductances (Gk). 


The resistance modulation of element k can be written : 
A^t = A T k = -^f A kl (A i + 2 RqjIqAI + Iq AR, j 


(Eq.ll) 



or as the matrix expression (with S kl = 1 for k = l , Su = 0 for **/.): 
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With the particular discretization illustrated in Fig.l, the heat dissipation in the segments 
can be written as 
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(Eq.14) 


A Prf,k ~ 


A 

2N 


k = 2,3, ...,N 


(Eq. 15) 


The matrix expression above provides the connection between the modulation in dissipated 
rf power and the modulation of device resistance. The first of the three terms in the 
expression for B kl represents the direct response to the //heating. The last term is due to 
the change in direct ohmic heating when the device resistance changes, while the middle 
term represents heating due to mismatch between the device and the intermediate frequency 
load resistor. For convenience, the following response parameter is defined: 



(Eq. 16) 


In the time domain, the //voltages over the device are 
V LO cos co LO t + V s cos co s t 

and cause a heat modulation at the //frequency of: 
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The resulting intermediate frequency voltage is: 
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(Eq. 18) 


(Eq. 19) 


and, finally, the single-sideband mixer conversion efficiency r\ can be formulated: 
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2.4 Thermal Fluctuation Noise 

A superconducting bolometer generates both Johnson (or thermal) noise and thermal 
fluctuation noise. This paper will only treat thermal fluctuation noise in any detail, although 
a more complete theory should also include Johnson noise. In practice, however, thermal 
fluctuation noise usually dominates for a diffusion-cooled superconducting HEB, when it 
is operated under the conditions that give minimum receiver noise. 

The origin of thermal fluctuation noise is the random exchange of heat carriers between the 
bolometer and the ambient temperature bath (the normal metal leads), and between different 
internal parts of the bolometer. In a system in thermal equilibrium with a large reservoir, 
thermodynamics gives the following well-known expression for the time variance of the 
energy of that system, where Cy is the heat capacity and kg is Boltzmann's constant: 

((A E{t)f) = k B C v T 2 (Eq.21) 

In a more relaxed nomenclature this relation can also be expressed as a fluctuation of the 
system temperature, although this is strictly speaking not correct since temperature is not a 
fluctuating parameter in a system in equilibrium: 

((A mf) = k B T 2 IC v (Eq.22) 

A real bolometric mixer is per definition not in equilibrium with the surrounding medium, 
since it relies on finite amounts of dc and rf heat being dissipated inside the device, and of 
that heat being transported out into the medium. However, if the inelastic electron-electron 
time is short enough, a local quasi-equilibrium can be established, where any small portion 
of the device only exchanges heat carriers with the parts of the device that are immediately 
adjacent to it. The exact conditions under which the expressions above no longer apply are 
hard to estimate, and in the end a mixer theory for thermal fluctuation noise will need to be 
compared to experimental results. 

For calculating the mixer noise, a circuit model is used, as shown in Fig.4, that is 
similar to the one used in the conversion efficiency calculation. The heat capacitances and 
heat conductances are the same as previously discussed, but the heat sources Qk are 
arranged differently. The noise heat sources represent exchange of energy between 
different segments, not dissipation, and are considered incoherent in time (white noise) and 
uncorrelated with each other. This is a Langevin approach that has similarities to some 
other formulations found in literature [19]. The method consists of two steps: first the 
magnitudes of the noise sources are determined, and then the amount of noise power 
coupled to the output load is calculated. 




Fig.4; The thermal circuit model used in calculating thermal fluctuation noise. 


The expression for the temperature modulation in the segments is (within a small 
frequency interval df): 
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The following formula is then used to calculate the source intensities, where the first 
equality follows from the Wiener-Khinchine theorem and the second one from 
thermodynamics: 
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Since the sources are time-incoherent it follows, again from the Wiener-Khinchine 
theorem, that their frequency-domain representation is frequency-independent. The source 
distribution can therefore be calculated as the solution to the equation system: 



(Eq.26) 


With the noise source distribution known, it is relatively straightforward to calculate the 
effect on the device resistance: 
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The power delivered to the output load in the small frequency interval ^becomes: 
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(Eq.29) 


3. Calculation Example 

As an example, a calculation of mixer conversion efficiency was made in order to 
demonstrate some of the features of our model. Although the method allows an actual 
experimental resistance-versus-temperature curve to be used, for simplicity the idealized 
curve in Fig.5 which contains features of actual RT curves was chosen for this example. 
Two superconducting transition temperatures can be seen in the figure (as is the case in 
measured RT curves), the higher of which occurs in the central main part of the bridge, 
while the lower one is present near the ends of the bridge. This is a simplified way of 
including the suppression of the critical temperature from the proximity effect, that occurs 
in the regions of the bridge that are contacted by the normal metal leads. The device 
dimensions used in the calculation are similar to those of our actual HEB mixers: Length 
162 nm, width 150 nm, thickness lOnm, and the heat conductivity 

G 0 / T 0 = 0. 1 (W/ mK 2 ) . The device length included the end regions of the device with the 
suppressed critical temperature (12 nm ). The ambient temperature was chosen to be 
4.24 K. 

Fig.6 shows the current- voltage (IV-) curves calculated for the device, with and without 
local oscillator power. The overall shape is similar to measured curves for our actual 
devices. At zero dc bias voltage, however, the curves become singular, since the model 
does not include the limitation on the dc supercurrent that is present in real superconductors 
[20], The mixer conversion efficiency shown in Fig.6 approaches a maximum at the point 
where the device cannot be biased in a stable way by the embedding circuit. This is also 
qualitatively consistent with typical experimental results. The circle in the figure marks 
approximately the best low receiver noise operating point that would be reachable in a real 
experiment. In addition, a peak in the conversion efficiency can be seen just above a bias 
voltage of 3 mV. This effect occurs when enough heat is dissipated in the device so that 
the end regions of the bridge heat up to their local (suppressed) critical temperature. This is 
sometimes observed as peaks in the coupled output noise in experimental data. Figure 7 
shows the calculated conversion effiency and output noise for the device when operated at 
the bias point marked in Fig.6. As seen, the conversion effiency curve has the expected 
single-pole shape, in this case with a 3 dB roll-off frequency of 3-4 GHz. The shape of 
the output noise curve, however, differs from that calculated from previous "lumped- 
element" theory [21], according to which it should have the same frequency dependence as 
that of the conversion efficiency. It should be pointed out that if the number of bolometer 
segments is reduced to one, the two theories become identical. Also, in a more realistic 
mixer model the Johnson noise should be included, which has not been done here. 




Fig.5; The R-versus-T curve used in the calculation is an idealized curve based on an 
experimental curve for a niobium HEB. 



Fig.6; Calculated IV-curves and calculated conversion effiency. 



Fig.7; Calculated conversion effiency and coupled output noise as a function of frequency. 
The noise here is defined as coupled power per Hz to the if load resistor, divided by kB- 





Summary 


A mixer model for diffusion-cooled hot-electron bolometers has been developed. A non- 
linear large-signal finite-difference method is used to calculate linearized device parameters, 
which are subsequently used in small-signal calculations of mixer conversion effiency and 
noise. The non-uniform temperature distribution inside the bolometer is taken into account 
both in the large- and small-signal analysis by discretizing the device into a large number of 
segments, and by assigning a local temperature to each segment. The small signal mixer 
analysis is done in the frequency domain, using matrix expressions to connect modulations 
in heating, temperature and resistance. Initial calculations using device parameters that are 
similar to those of real bolometers give values for conversion and noise that seem 
reasonable. However, a detailed comparison to experimental data still remains to be done, 
and will be presented in a future publication. Future work will also include a model for 
Johnson noise. 
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